Prospects for direct detection of circular polarization of gravitational-wave 

background 
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We discussed prospects for directly detecting circular polarization signal of gravitational wave 
background. We found it is generally difficult to probe the monopole mode of the signal due to 
broad directivity of gravitational wave detectors. But the dipole (I = 1) and octupole (I = 3) modes 
of the signal can be measured in a simple manner by combining outputs of two unaligned detectors, 
and we can dig them deeply under confusion and detector noises. Around / ~ O.fmHz LISA will 
provide ideal data streams to detect these anisotropic components whose magnitudes are as small 
as ~ 1 percent of the detector noise level in terms of the non-dimensional energy density flGw(f)- 

INTRODUCTION 



As gravitational interaction is very weak, significant efforts have been made to detect gravitational waves. But, 
on the other hand, we will be able to get rich information of the universe by observing gravitational waves that 
directly propagate to us with almost no absorption. Various astrophysical and cosmological models predict existence 
of stochastic gravitational wave background, and it is an interesting target for gravitational wave astronomy Q] . For 
its observational prospects, we need to understand how we characterize the background and what aspects we can 
uncover with current and future observational facilities. 

One of such aspects is circular polarization that describes whether the background has asymmetry with respect 
to magnitudes of right- and left-handed waves. Circular polarization of gravitational wave background might be 
generated by helical turbulent motions (see e.g. Q and references therein). Inflation scenario predicts gravitational 
wave background from quantum fluctuations during acceleration phase in the early universe, but asymmetry of 
0^ ' left and right-handed waves can be produced with the gravitational Chern-Simon term that might be derived from 
string theory and might be related to creation of baryon number (see e.g. 0] and references therein). Primordial 
gravitational wave background, including information of its circular polarization Q, can be indirectly studied with 
CMB measurement pj at very low frequency regime / ~ 10 _17 Hz that is largely different from the regime directly 
' accessible with gravitational wave detectors studied in this Letter. Gravitational wave background from Galactic 
binaries can be polarized, if orientation of their angular momentum have coherent distribution, such as, correlation 
with Galactic structure. While observational samples of local Galactic binaries do not favor such correlation [(| , this 



will also be an interesting target that can be directly studied with LISA. 



It is well known that observation of gravitational wave is intrinsically sensitive to its polarization state [2j . This is 
because we measure spatial expansion and contraction due to the wave, and polarization determines the direction of 
the oscillation perpendicular to its propagation direction. Another important nature of the observation is that we have 
to simultaneously deal with waves basically coming from all the directions, in contrast to observations with typical 
electromagnetic wave telescopes that have sharp directivity. Therefore, polarization information and directional 
information couple strongly in observational analysis of gravitational wave background. In general, orientation of 
a gravitational wave detector changes with time, and induced modulation of data stream is useful to probe the 
polarization and directional information. For ground based detectors, such as, LIGO, this is due to the daily rotation 
of the Earth ||. For a space mission like LISA, this change is determined by its orbital choice jj|E3- In addition, we 
can also expect that several independent data streams of gravitational waves will be taken at the same time [ill Il2| . 
In this Letter, in view of these observational characters, we study how well we can extract information of circular 
polarization of the background in a model independent manner about its origin. 



FORMULATION 



The standard plane wave expansion of metric perturbation by gravitational waves is given as 

h ab (t,x)= r df I ^M/,«)e 2 ™ /( '- n:E) ef b (n), (1) 

where S 2 is the unit sphere for the angular integral, the unite vector n = (sin 9 cos 0, sin 8 sin <f>, cos 9) is propagation 
direction, and e+ fc and e* b are the basis for the transverse-traceless tensor. We fix them as e^ b = eg (g> eg — CS> and 
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e* b = eg(g)e. < f 1 + e ( f ) ®eg where eg and are two unit vectors with a fixed spherical coordinate system. As the metric 
perturbation h a b(t, x) is real, we have a relation for complex conjugate; hp(—f, n) — hp(f, n)*. When we replace the 
direction n — > n, the matrices have correspondences e+ h — > e+ h (even parity) and e* b — > — e* h (odd parity). 

To begin with, we study gravitational wave modes at a frequency /. Here, we omit explicit frequency dependence 

for notational simplicity, unless we need to keep it. The covariance matrix ( fcWSWl) ^ 

two polarization modes h + (n) and h x (n) is decomposed as 

5 drc (n-n r ) ( I + Q U-iV\ ^ 



hUn)h x (n')) (h x (n)h* x (n>) 
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where the symbol (. . .) represents to take an ensemble average for superposition of stationary incoherent waves, and 
5drc{ ) is the delta function on the unit sphere S 2 . We have defined the following Stokes parameters 01 analog 
to electromagnetic waves as I(n) = (\h + \ 2 + \h x | 2 ) /2, Q{n) = (\h + \ 2 - \h x | 2 ) /2, U(n) = (h + h* x + h%h % ) /2, 
and V(n) = i(h+h* x — /i5_/i x ) /2. The parameter I(n) represents the total intensity of the wave, while Q(n) and 
U(n) is related to linear polarization. The parameter V(n) is related to circular polarization, and its sign shows 
whether right- or left-handed waves dominate. When we rotate the basis vectors eg and around the axis n 
by angle ip, the parameters / and V are invariant (spin 0). But the parameters Q and U are transformed as 
(Q ± iU)'(n) = e^ 4 ^(Q ± iU)(n), and the combinations (Q ± iU)(n) have spin ±4 [l]f. 

Next we discuss response of a two-arm interferometer J with 90° vertex angle and equal arm- length L. We put l± 
and I2 as unit vectors for the directions of the arms. The beam pattern functions Fj = (l\ ■ e p ■ l\ — I2 ■ e p ■ l2j/2 
(P = +, x) represent relative sensitivities to two linearly polarized gravitational waves (P = +, x) with various 
directions n 0- Note that the function Fj~ has even parity and Fj has odd parity with respect to the direction n 
(see e.g. Q). 

In this Letter we mainly deal with low frequency regime with ///* -C 1. Here, using the arm-length L, we have 
defined a characteristic frequency /* = l/(2wL) that corresponds to lOmHz for LISA and 1Hz for Big Bang Observer 
(BBO) [HG]!. LISA is formed by three spacecrafts that nearly keep a regular triangle configuration 9]. From 
its six one-way data streams, we can make Time-Delay-Interfcrometer (TDI) variables that cancel laser frequency 
noises. We can select three TDI variables A, E and T whose detector noises are not correlated and can be regarded 
as independent At the low frequency regime the responses of the A and E modes can be effectively regarded 

as those of two-arm interferometers whose configuration are shown in figure 1 In this figure we put the whole 

system on AF-plane (8 = tt/2). The beam pattern functions of the A mode are explicitly given as 



F^(0,<j)) = -(1 + cos 2 61) cos 20, F%(8,4>) = - cos sin 20, (3) 



while those for the E mode are given by replacing <f> — > (f> — ir/4 in the above expressions. The beam pattern functions 

p -m 4 )(i 
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for the T mode are quite different from the A and E modes, and given as F^ = J^ i=1 (m; • e F ■ m,)(mi • n) (P = +, x ) 
where directions of three unit vectors are shown in figure 1 |l4t Il7| . At low frequency regime sensitivity of the T 
mode to gravitational waves is ~ (///*) _1 times worse than those for the A and E modes 0,^3- Therefore we put 
this mode aside for a while. 

We can express responses of the A and E modes to gravitational waves from a single direction n as ( J = A, E) 



rj(n) = (F+h + (n)+F*h x (n)) [l + i{f/U)Dj{n) + 0(f 2 /f 2 )] , (4) 

where the second factor [. . .] is corrections caused by the finiteness of the arm-length, and two real functions Da(&, <j>) 
and De(0,4>) depend on the propagation directions n of waves (see e.g. ^3)- I* 1 eq.JU we neglected overall factors 
that depend only on frequency / and are irrelevant for our study. These formal expressions for the perturbative 
expansion with respect to the ratio (///*) can be generally used with relevant beam pattern functions, including the 
case for the T-mode or Fabri-Perot detectors as LIGO. 

The information V cannot be produced from the response or alone. To get it we need independent linear 
combination of h + and h x . For example the term (r a{0 , 4>) r a(9 , 4 1 )*) can be written only with I,Q and U. We now 
take the low frequency limit, and keep only the leading order terms of the expansions. The cross term (rA{n)rE{n')*) 
is given as 
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FIG. 1: The effective two-arm interferometers corresponding to the TDI modes A and E. Three spacecrafts of LISA are shown 
with circles that are put on the XY-plane (9 = n/2). The beam pattern functions for T mode are given with three unit vectors 

rra as = J2^=i( mi ' eP ' m i)i m i • n ) (P = +, x). 
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With the following combination 
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we can extract the circular parameter V alone. As shown in figure 1, the effective detector E is obtained by rotating 
the detector A around Z-axis by tt/4 [13. If this angle is 6, there appears a factor sin 26 in the final expression in 
eq.©. Therefore, in some sense, LISA will provide an optimal set (A,E) to study the parameter V. In contrast, 
sensitivity of correlation analysis to the monopolc intensity I 00 is proportional to cos 2(5, and LISA cannot probe it 
with the method, as is well known. 

We now discuss responses of detectors to gravitational waves from all directions n. Firstly, what we can get obser- 
vationally form the A and E modes are the following integrals Ra — J rA(n)dn and Ra — J rA(n)dn. Considering 
the spin of quantities /, Q HU and V we can expand them in terms of the spin weighted spherical harmonics s Yi m as 
I(9,<f>) = J2limYi m (9^), (Q±iU)(9,4>) = EK±,irn-±4Y lm (9,ci)), V{9,4>) =EVi m Yi m (9,<f>), where the coefficients 
K±,im are defined for I > 4. With these harmonic expansions and eq.©, the combination C — Im (RaR*e) is evaluated 
as C = — 



becomes 



yV-jo. When we rotate two interferometers A and E with Euler angles (a,/3, 7), the signal C 



C(a,/?,7) 



16tt 
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V lm Y lm (a, (3) - —V 3m Y 3m (a,(3). 



(7) 



This result is obtained by relating the coefficients Vi m in original coordinate with those in rotated coordinate |8|]. 
LISA moves around the Sun with changing the orientation of its detector plane that is inclined to the ecliptic plane by 
60 degree 0- To describe these motions in the ecliptic coordinate, the Euler angles are given as a — 27r(t/lyr) + a Q , 
(3 = 7r/3 and 7 = — 27r(t/lyr) + 70 with time t and constants ao and 70 ^H]. This parameterization is also valid 
for BBO ^^|. The parameter 7 determines the so-called cartwheel motion, but the combination C(a,(3, 7) does not 
depend on it. From eq.(0 we can understand that the observed signal C(a, {3, 7) with LISA is decomposed into 
modulation patterns with frequencies / = 0, ±1/3, ±1/2 and ±lyr _1 from m = 0, ±1, ±2 and ±3 modes respectively. 

So far we have used the low frequency approximation with ///* <C 1. When we increase the frequency /, the 
correction terms i(f / f*)DA(0,<j)) and i(f / f*)DE(0,(f>) in eq.Q change the phases of rA{9,<j)) and rE(9,<fr) as a 
function of propagation directions n. With these correction terms, the combination RaR*e ~ R\Re has contributions 
of /, Q and U modes, and we cannot extract the circular polarization V in a clean manner. Note also that the first 
order term O (///*) for the combination RaR* e + R\Re depends on the parameter V due to the corrections. We can 
restate the situation as follows; Roughly speaking, the circular polarization is measured by correlating two data with 
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phase difference ir/2 0. This is related to the fact V(n) oc Re {ta(^, <j))rE(d, </>)*e i7r / 2 ). But the finiteness of the 
arm-length modulates the phase as a function of direction, and we cannot keep the phase difference 7r/2 simultaneously 
for all the directions. When we consider the spatial separation d of two interferometers, same kind of arguments hold 
by perturbatively expanding the phase difference exp[27ri/e£ • n] with a expansion parameter in addition to 

that with (///*) for the effects of arm-length. Therefore, in our analysis, the requirement for the low frequency 
regime is not for simplicity of calculation, but is a crucial condition to extract the circular polarization alone with 
gravitational wave detectors that have broad directional response and are remarkably different from electromagnetic 
wave telescopes. 

As shown in eq.©, we cannot measure the monopole moment Vbo of circular polarization using the signal C made 
from the A and E modes due to parity reason at low frequency limit. More specifically, products, such as F^F^ , 
have odd parity. This simple results hold for signal C with any two two-arm interferometers at low frequency limit. It 
does not matter whether their vertex angles are not ir/2, whether two arm-lengths of each detector are not equal, or 
whether two interferometers are on unparallel planes (e.g. for LIGO and VIRGO combination). The signals C with 
(A, T) or (E, T) modes depend on the coefficients Vi m with even / at their leading order. But, in the case of LISA, the 
combinations do not have monopole moment Voo- This is because of the apparent symmetry of these data streams 
[l2l Hij]. Furthermore we cannot get the moment Voo even using their higher order terms with 0((///*) n ) (n > 1), 
as long as LISA is symmetric at each vertex. A future mission might use multiple LISA-type sets with data streams 
{A, E,T},{A', E',T'}, ... ^J. We confirmed that even if detector planes for A and T" modes are not parallel, the 
monopole Voo can not be captured by the signal C with their combination at their leading order. 

As we discussed so far, it is not straightforward to capture the monopole Voo m a simple manner. But, in principle, 
we can manage this. For example, we add a detector E% that is given by moving the original detectors E in figure 
1 by distance d toward +Z-direction, and then take the signal RaR*e 2 - At order 0((/rf) 1 (///*) ), we can probe the 
monopole Voo (3- This kind of arrangement for detector configuration might be important to study anisotropy and 
polarization of gravitational wave background in the long run. 



OBSERVATION WITH LISA 



Next we discuss how well we can analyze the information V; m of circular polarization of stochastic gravitational 
wave background with LISA. The observed inclinations of local Galactic binaries are known to be consistent with 
having random distribution with no correlation to global Galactic structure Q. This suggests that the confusion 
background by galactic binaries is not polarized. In future, LISA itself will provide basic parameters for thousands of 
Galactic binaries at / > 3mHz, including information of their orientations [9j. This can be used to further constrain 
the allowed polarization degree of the Galactic confusion noise. While this information will help us to discriminate 
the origin of detected polarization signal with LISA, our analysis below does not depend on this outlook. Our target 
here is a combination pQcw,b(f) where Qcw,b(f) is the normalized energy density of gravitational wave background 
(see for its definition), regardless of its origin. The parameter p = (\/J2 m l^iml 2 + S m \^m\ 2 )/It is its circular 
polarization power in I = 1 and 3 multipoles (It = \/X^m \^im\ 2 '■ total intensity). For simplicity, we do not discuss 
technical aspects in relation to the time modulation of the signal C due to the motion of LISA that was explained 
earlier around eq.ljZJ). We can easily make an appropriate extension to deal with it [Toj ]. 

As we will see below, the signal C is a powerful probe to study the target p&>GW,b(f) in a frequency regime where 
the amplitude p£law,b(f) itself may be dominated by the confusion background noise or detector noise. We take 
a summation of the signal C for Fourier modes around a frequency / in a bandwidth A/ ~ /. As the frequency 
resolution is inverse of the observational time T t s and the total number of Fourier modes in the band is (T Q b s Af), 
the expectation value for the summation becomes C(T f, s A/) ~ pI t (T b s Af). Here we used the fact that the detector 
noises for the A and E modes are not correlated. In contrast the fluctuation N for this summation is given as 
N = (riA(f)nE(f)* — n A(f)*nE(f)) , where the total noises n^(/) and Ue(J) include both detector noise and 

1/2 

circularly unpolarized potion of the confusion noise. The root-mean-square value of N is written as ^| A^| 2 ^> ~ 
S n (f)(T i JS Af) 1 ^ 2 with S n (f) being the total noise spectrum for the A and E modes. Thus the signal to noise ratio for 
the measurement is SNR ~ plT(T obs Aff/ 2 / S n (f). We can rewrite this expression with using the normalized energy 

density flow and obtain SNR ~ ( jncw j (TobsAf) 1 / 2 , where the magnitude flcw.n corresponds to the total noise 
level around frequency /. When the detector noise is dominated by the background noise, we have £low,n ~ &GW,b- 
This will be the case for LISA around 0.3mHz < / < 2mHz. As a concrete example, we put / ~ O.lmHz with the 

expected noise level fJ G w> ~ 1CT 10 9], and take bandwidth A/ - /. Then we have SNR - ( ^"j" ) (T/Syi) 1 / 2 . 
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The improvement factor (T (, S A/) x / 2 for the detectable level of the target pQcw,b is caused by the same reason as 
standard correlation technique for detecting stochastic background 0. Interestingly, this factor (T obs A/)" 1 / 2 - 0.01 
around / ~ O.lmHz is almost same as the maximum level of relative contamination (/ If* ) ~ 0.01 for LISA by other 
modes (I, Q, U) due to finiteness of arm-length. A Japanese future project DECIGO [20j is planed to use Fabri-Perot 
type design with its characteristic sensitivity /* = 50Hz [2l|, while its best sensitivity is around ~ 0.3Hz similar to 
BBO whose characteristic frequency is /» = 1Hz. Therefore, DECIGO is expected to be less affected by other modes 
and has potential to reach the level pQcw.b ~ 10~ 16 with one year observation. 

If circular polarization parameter V(0, <j>) of cosmological background is highly isotropic in the CMB-rest frame, 
its dipole pattern will be induced by our peculiar motion [2^. This might be an interesting observational target in 
future. In this case we can predict the time modulation pattern C(a,P,j) for rotating detectors. A consistency check 
using this modulation will be powerful approach to discriminate its cosmological origin. Actually a similar method 
can be used for LISA to analyze circular polarization of Galactic confusion noise. 
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